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contraction and variational techniques with the real-space renormalization 



group approach. It applies to lattice systems of infinite extent and is ideal for 



studying phase structure and critical phenomena. The CORE approximation 
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Many problems in particle and condensed matter physics cannot be studied with conven- 
tional perturbation theory. Aside from Monte Carlo simulations, few tools allow one to deal 
with general Hamiltonian systems and fewer tools deal directly with their infinite- volume be- 
havior. This paper introduces a new tool, the Contractor REnormalization group (CORE) 
approximation, which can handle this class of problems. The CORE approach is a simple, 
systematic procedure for improving any Hamiltonian real-space renormalization group cal- 
culation. Its virtues are: it is a variational procedure which is systematically improvable to 
any desired degree of accuracy; it applies to lattice systems of infinite extent, allowing direct 
study of phase structure and critical phenomena; it provides tools for error estimation; it 
requires modest computer resources by modern standards; it is complementary to Monte 
Carlo methods; systems with dynamical fermions can be treated. 

We start with a brief description of the CORE approximation, then illustrate and test 
the method in two different applications to the 1+1-dimensional Ising model. 

Basic Ideas Choosing a good trial state is crucial to the success of any variational 
calculation, especially one involving a large number of degrees of freedom. The Hamiltonian 
real-space renormalization group (RSRG) method [JTJ] is an algorithm for constructing a 
class of trial states appropriate for lattice systems. In this approach, one partitions the 
lattice into blocks containing a few sites and diagonalizes the Hamiltonian associated with 
each block. One then thins the Hilbert space by discarding all states except those which 
are tensor products of some chosen subset of low-lying block eigenstates, and an effective 
Hamiltonian which describes the mixing of the retained states is computed. This truncation 
process is iterated until the effective Hamiltonian evolves into a fixed form which can be 
easily diagonalized. 

Unfortunately, simple RSRG truncation procedures tend to severely underestimate the ef- 
fects of block-to-block couplings and this hinders the accurate description of long- wavelength 
modes on the full lattice. Past approaches to overcoming this problem have concentrated on 
using larger blocks, increasing the number of states retained per block, or introducing more 
sophisticated truncation schemes @,[|. The t-expansion has also been used f§]. The CORE 

2 



approximation is a new approach to this problem which emphasizes simplicity, versatility, 
and insensitivity to the precise details of the truncation scheme. This insensitivity, or ro- 
bustness, frees one from the need to develop clever truncation algorithms or to retain many 
states per block. This feature of CORE, combined with its simplicity, greatly enhances 
its usefulness for higher dimensional systems. CORE also allows the use of manifestly 
gauge-invariant RSRG schemes when studying lattice gauge theories; such simple schemes 
cannot be exploited in the naive multi-state approach since gauge non-invariant states are 
necessary for coupling neighboring blocks after the first truncation step. Furthermore, the 
CORE method does not suffer from the series reconstruction difficulties which plague the 
t-expansion. 

The basic idea of the CORE approach is to use contraction techniques to steer the RSRG 
iteration. In the limit t — > oo, the operator e~ tH contracts any trial state |$ V ar) onto the 
lowest eigenstate of H with which it has a non- vanishing overlap. Therefore, the expectation 
value 
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tends to the corresponding eigenvalue eo of H as t becomes large. In general, £{t) cannot be 
computed exactly. Reliably approximating £{t) is an integral part of the steering process in 
the CORE method. 

An important step in building a CORE approximation to £(t) is to construct an easily 
computable operator T(t) which closely approximates e~ tH for t in some range < t < t max . 
To find such an operator ||, first divide H into two (or more) parts, i.e., H = Hi + H2, 
where the individual parts Hi and H2 are chosen such that e~ tHl and e~ tH2 can be computed 
exactly. Next, rewrite e~ tH as a symmetric product 

e -tH _ g-tffi/2 e -tH 2 /2 e C 3 (t) e -tH 2 /2 g-tETi/2^ ( 2 ) 

where C${t) is a sum of terms all of which begin in order t 3 or higher. The simplest T(t) 
is obtained by replacing e° 3 ^ by the identity operator. One way to construct a better 
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approximation is to retain low-order terms in Cs(t) and rewrite the exponential of these 
operators as a symmetric product of explicitly computable terms. Another is to use the 
operators T p (t) = [T(t/p)] p . In any case, it is very important to ensure the approximate 
contractor satisfies all the symmetries of H. 

Given a contractor T(t), e can then be bounded from above by computing 
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A best estimate for eo is obtained by minimizing £t(£) with respect to t and any parameters 
in |$ var ). For a trial state |$ var ) = YJj=\ a j\4>j)i where is some set of orthonormal 

states, one can show that minimizing £r(i) with respect to the a.j parameters is equivalent 
to solving the generalized eigenvalue problem 

det([T(t)i/T(t)]-A[T(t) 2 ]) =0, (4) 

where [. . .] denotes truncation to the subspace spanned by the \<pj) states. In particular, for 
an operator O, [O] = POP^ where P is the projection operator P = £)?=i Thus, 
finding the best trial state |$ V ar) is equivalent to diagonalizing the effective Hamiltonian 

H cS (t) = [TitfY 112 {T{t)HT{t)\ [T(t) 2 ]- 1/2 - (5) 

Developing this operator in the RSRG iteration instead of [ H ] is a key innovation of the 
CORE approach. 

The effective Hamiltonian defined by Eq. [| cannot be exactly computed. Another novel 
feature of the CORE approach is the use of the finite cluster method to evaluate H e ^{t). 
In this method, H e s(t) (or any other extensive quantity) is calculated as a sum of finite- 
volume contributions (see Ref. 0). The finite cluster method, which will be described 
later when applying CORE to the Ising model, is simple to implement, provides numerous 
computational checks, and does little or no harm to the variational bound in Exit)- 

The final ingredient in the CORE approximation is the selection of a best value for t in 
each RSRG step. This can be done in a number of ways. One can extract the coefficient 
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of the identity operator in H e g and vary t to minimize this quantity. Better yet, one can 
evaluate H eS in a simple product state to produce a mean-field estimate of the ground state 
energy and minimize this with respect to t. 

In summary, the CORE method generates a sequence of effective Hamiltonians ii^g (i* ) 
by successive thinning of degrees of freedom using the recursion relation 

H^ +1 \t) = Rn(t)lT {n \t)H^\r n )T^(t) ]R n (t) } (6) 

where R n (t) = [ T^ n \t) 2 J^ 1 ^ 2 and the contractor T^ n \t) is constructed to approximate 
exp[— tH$(t*)]. Eq. | is evaluated using the finite cluster method and a best t — t^+i 
must be chosen. As the recursion proceeds, the effective Hamiltonian evolves eventually 
into a simple form which can be trivially diagonalized, yielding estimates of the ground 
state energy and the energies of some low-lying excited states. 

The expectation value of an extensive operator O can also be evaluated in the CORE 
method. One develops O using the same RSRG transformations as for H, producing a 
sequence of effective operators 0^(t*). The matrix element of O e ff is then evaluated once 
H e g has evolved to the point where its ground state can be easily determined. 

Note that from a programming point of view, CORE calculations involve mainly matrix 
multiplications; diagonalizations and inversions of only very small matrices are required. 
Often the matrices will be sparse and one can exploit efficient algorithms for multiplying 
them. 

The 1 + 1-Dimensional Ising Model We illustrate and test the CORE approximation in 
two different applications to the 1+1- dimensional Ising model. The Hamiltonian in this 
model is given by 

Rising = - J2 [ c a°" z (j) + Sxa x (J)a x (j + 1)] , (7) 
j 

where j labels the sites in the infinitely long chain, c\ = cos(A7r/2), and s\ = sin(A7r/2), for 
< A < 1. This model exhibits a second-order phase transition at A = 1/2. For A < 1/2, the 
ground state is unique and the order parameter (a x (j)) = 0, for some site j. When A> 1/2, 
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the ground state is twofold-degenerate and the order parameter takes values (<r x (j)) = 
±(1 -cot 2 (A 7 r/2)) 1 /8. 

The CORE approximation is best applied in the following sequence of steps: (1) choose 
an RSRG algorithm by specifying how to partition the lattice into blocks and which states 
to retain on each block; (2) specify the truncation order in the cluster expansion of H c{i ; (3) 
deduce the general form of H e g based on the choices made in steps 1 and 2; (4) construct a 
contractor T(t) which closely approximates exp(— tH e s) and is easily computable; (5) choose 
a method of determining the optimal value of t in each RSRG step; (6) iteratively compute 
H e s using Eq. || with initial condition = H, where H is the Hamiltonian of interest, 
until H e{i can be easily diagonalized. 

In our first application, we partition the lattice into two-site blocks and truncate the 
Hilbert space to the lowest two eigenstates in each block. Since our intention here is to 
carry out only the simplest of calculations, we choose to truncate the cluster expansion 
of H e fi(t) after three-block clusters. The general form of H eS (t) may then be deduced by 
considering how it is computed in the finite cluster method. 

Evaluation of H eS (t) by the finite cluster method is accomplished in the following se- 
quence of steps. First, compute H e{i (t) using Eq. ^| on a sub-lattice which contains only a 
single block. This yields 

h { S(t)=H$(t) = cV(t)u + c^(t)a zl (8) 

where u is a 2x2 identity matrix. Next, calculate H e g(t) for a theory defined on a sub-lattice 
made up of two adjacent blocks. This Hamiltonian takes the form 

H%® = c^\t)u L u R + cf\t){a L z u R +u L a*) 

+ <& (* )°M+<S (*)°M+C$ (t)aX, (9) 

where L and R refer to the left and right blocks, respectively, in the cluster. Remove from 
the two-block calculation those contributions which arise from terms already included in the 
single-block calculation: 
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h%(t) = H^{t)-u L ® h^ R \t)-h^ L \t)®u R . (10) 

Repeat this procedure for sub-lattices containing successively more connected blocks, then 
sum the contributions h^\t) from these sub-lattices with weights given by the number of 
ways each sub-lattice can be embedded in the full lattice. The stage at which one cuts off 
this cluster expansion determines the maximum range of the interactions which will appear 
in H eS . For our choices, the effective Hamiltonian in this model takes the general form 

H cS (t) = -£c Q (t)O a (j), (11) 
O a {j) = a ao (j)a ai (j + 1) . . . a ar (j + r), (12) 

where c a (t) are the couplings, a labels the different types of operators, and j is a site label. 
There are only two one-site operators: = {u, z}, where u denotes the identity operator. 
In other words, the only one-site operators are O u (j) = a u (j) = 1 and O z (j) = <7 z (j)- 
There are three two-site operators: = {xx,yy, zz). The three-site operators are cS^ = 
{xzx, xux, xxz, zxx, yzy, yuy, yyz, zyy, zuz, zzz}. 

Our first contractor is built using the approximation exp [— tH e g(t*)] ~ S^fjSft), where 
the operator S(t) = Y\ a {Y\j exp[tc a (t*)O a (j) /2]}. The operators in the a product are ordered 
according to their site range, increasing in size from right to left. This operator can be 
simplified using exp[yO a (j)] = cosh^y [1 + tanhy O a (j)}, where iV is the (infinite) number 
of sites in the lattice. Discarding the unimportant cosh 7 ^ factors, one obtains a contractor 
given by T x (t) = S^S^i), where S 1 (t) = UaiUj [1 + tanh(c a t/2)O a (j)}}. 

Lastly, t is chosen in each RSRG step to minimize the expectation value of H e g(t) 
evaluated in the mean-field state given by \ip m f) = lT,( cos $| Tj) + e^sin#| Lj)), where 
a z(j)\ Tj) = | Tj) and (T z {j)\[j) = — | The matrix element (V'mfl-Heff(OlV'mf) is minimized 
with respect to t, 9, and 0, simultaneously. 

For our second application of the CORE method, the lattice is divided into blocks con- 
taining three sites and the Hilbert space is again truncated to the lowest two eigenstates 
in each block. The cluster expansion is taken only to three-block clusters, so H ef[ takes 
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the general form shown in Eqs. |Tl| and [12|. We use an approximate contractor given by 
T 2 (t) = S , 2(t)S , 2 (t) with S 2 {t) = exp(—tV/2)exp(—tH b /2), where H b contains all intra-block 
interactions and V contains all inter-block operators (those which cross block boundaries). 
Note that exp(-tH b /2) = U p exp(-tH b (p)/2) and exp(-tV/2) = U p exp(-tV(p)/2), where 
H b (p) contains all operators which solely act on block p and V(p) contains only interactions 
between block p and p+1. The operators H b (p) and V(p) can be exponentiated numerically 
with no difficulty We fix t by minimizing the expectation value of H e g in the mean-field 
state as described previously. 

Selected estimates Eq of the ground-state energy density from both variants of the CORE 
approach described above are compared to the exact energy density e in Fig. [L| Calcu- 
lations were done using T™{t/ri) and T£(t/n) for various values of n. The fractional errors 
5 eo shown in this figure are defined by 5 eo = \(E — eo)Ao|- Selected mass gap estimates A 
are compared to the exactly-known gap in Fig. 0. Fig. |3| illustrates the amounts by which 
the T 2 12 CORE estimates of the magnetization M. = \(a x (j))\, for some site j, differ from 
the exact values. The accuracy of the results is striking, especially considering that only 
the first three terms in the cluster expansion were included in the calculations. The CORE 
method reproduces the correct location of the critical point with remarkable precision. The 
critical exponent ( was extracted from a straight-line fit of our T 2 12 results for ln.M to the 
form ln.M = C l n (l — Ac/^ 2 )' where A = tan(A7r/2), A c = tan(A c 7r/2), and A c is our com- 
puted value for the critical point. For A c = 0.5053 and fitting in the range 0.51 < A < 1.0, 
we obtain ( = 0.12437, to be compared to the exact value of 0.125. The CORE procedure 
produces better results for a given effort than multi-state RSRG methods previously used 
||. The results also compare very favorably to previous t-expansion calculations 0. Using 
larger blocks or including more terms in the cluster expansion should further improve these 
results. 

Conclusion Given its simple theoretical foundations, the relative ease of implementa- 
tion, and our success in applying it to the 1+1-dimensional Ising model, we believe that the 
CORE approximation will prove to be a powerful tool for analyzing intrinsically nonpertur- 
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bative systems. One particularly exciting feature of this method is that it can be applied to 
systems containing dynamical fermions, systems which resist treatment by present stochas- 
tic means. In general, we feel that the possibility of eliminating the quenched approximation 
in lattice quantum chromodynamics, better studying spontaneous symmetry breaking and 
other nonperturbative phenomena in relativistic field theories, and probing the low-energy 
physics of the Hubbard and t — J models warrants further work with the CORE approxi- 
mation. 

This work was supported by the U. S. DOE, Contract No. DE-AC03-76SF00515, NSERC 
of Canada, and the UK SERC, grant GR/J 21347. 
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[G. 1. Fractional error <5 eo in the ground-state energy density estimates against A. Results 
T-j 2 (dashed curve), T x 16 (solid), and T 2 12 (diamonds with dotted curve, to guide the eye) are 
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FIG. 2. Mass gap estimates A against A. The diamonds and squares indicate CORE estimates 
obtained using 

T 16 

and T 2 12 , respectively. The exact mass gap appears as a solid curve. 
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FIG. 3. Magnetization M. against A. The diamonds indicate CORE estimates obtained using 
T 2 12 , the solid curve shows the exact magnetization, and the dot-dashed curve shows the estimates 
from mean-field theory. 



14 



